Kinetic study of hydrogen sulfide decomposition on Pt(111) surface

In this work, kinetic of H₂S conversion to H₂ molecule on the surface of Pt(111) is studied using kinetic Monte Carlo simulation. The results of simulation were fitted to the experimental temperature-programed desorption spectra. The good agreement between the empirical and the simulated data confirms the proposed mechanism and kinetic data (activated energies and pre-exponential factors). The influence of variables such as temperature and concentrations of H₂S and H₂ on the overall results of hydrogen production is studied. The condition is proposed in which the best yield of reaction at minimum temperature is obtained. Results show that platinum is a perfect catalyst for converting H₂S to H₂ and it has a perfect performance (98%) after 5 μs at low temperature of 227°C.     

Experimental and molecular simulating study on promoting electrolyte-immersed mechanical properties of cellulose/lignin separator for lithium-ion battery

Lithium-ion batteries have been developing intensively and earn an unprecedented reputation, yet advanced performance and safety issue still require considerable investigation. Separator is vital to comprehensive properties of batteries, where the mechanical properties are key to breaking through of new-type separator. Unfortunately, electrolyte submersion has caused damage to strength of cellulose separator. Whereupon, in this work, cellulose separator is optimized by introducing lignin particles to promote electrolyte-immersed mechanical strength. Experiments are conducted concerning surface morphology, contact angle, porosity, electrolyte uptake, mechanical properties and electrochemical performance. Molecular simulation is implemented to explore the mechanism of tensile behavior of cellulose and lignin subjected to electrolyte solvents. Experimental results confirm positive effect of lignin addition in improving mechanical properties and simultaneously maintaining impressive electrochemical performance of the cellulose/lignin composites separators. Besides, lignin addition amount of 2.5% and 5% is recommended to achieve promising overall properties. Molecular simulation has successfully unveiled that weakening of cellulose separator submerged in electrolyte is resulted by the deformed cellulose amorphous region and the promoting effect of adding lignin is contributed from the new hydrogen bonds generated between cellulose and lignin molecules. Hopefully, this work provides novel insight on preparing remarkable separator and mechanism of materials behavior.     

秋千参数自激振动的数学模型与数值模拟

将秋千视为变长度单摆,基于归一化高斯函数构造了摆球相对摆杆的速度表达式.在此基础上建立了秋千的非线性动力学微分方程,并据此数学模型对秋千运动进行了数值模拟,得到了秋千的运动规律和功能转化的机制.     

含水合物沉积物多相渗流特性研究进展

天然气水合物作为一种储量大、无污染的清洁能源近些年受到了广泛关注. 近20年来,中国进行了较大范围的陆海域天然气水合物储层勘探与储量预测.2017年,中国地质调查局牵头对南海神狐海域的天然气水合物进行了基于降压渗流原理的试验性开采.国内外已进行的水合物试采工程面临着气体产量低、出砂较多等问题,其最主要的原因之一是开发过程中沉积物内复杂多相渗流机理尚不明晰.本文综述了平行毛细管模型、Kozeny模型等广泛应用于天然气水合物开发渗流分析的理论模型,对比分析了水合物开发多尺度渗流过程模拟方法,简述了国内外含水合物沉积物渗透率测试、渗流过程中沉积物物性演变以及水合物开采室内模拟等方面的渗流实验进展,总结了矿场尺度的天然气水合物储层开采过程中产气数值模拟手段,展望了多相渗流模型、储层原位含水合物样品室内测试及结构与物性演化、矿场尺度数值模拟与水平井压裂技术等应用研究的未来方向与挑战.      

沙丘背风侧不同铁路结构形式对风沙环境的适应性分析

当铁路穿越大风沙漠地区时,风沙灾害会对铁路工程及其正常运营产生严重威胁,而设计一种合理的铁路结构形式能够减小风沙沉积对铁路工程的危害. 在本文中,以敦煌至格尔木铁路沙山沟段落为研究对象,采用多相流的方法对越过沙丘的风沙运动过程进行数值模拟,分别讨论了风沙运动对位于沙丘背风坡的铁路路基工程和桥梁工程的影响. 主要的模拟结果显示:路基工程明显降低了风速并且将沙丘后的回流区分成了两部分,而桥梁工程的导流效应则压缩了沙丘背风坡的回流区;轨道间的道碴增大了铁路表面的粗糙度,在轨道间有少量沙粒沉积,而路基工程两侧则有大量积沙;铁路表面的积沙量与摩阻风速呈现出非线性关系,随着摩阻风速的增大,路基工程沙粒沉积的增加速度大于风蚀能力的增加速度,而桥梁工程则正好相反. 在防止风沙危害铁路方面,设置桥梁工程明显优于路基工程. 本研究为风沙运动对铁路工程的影响提供了理论支持,也为今后的铁路工程设计提供了新的思路与研究工具.      

HyperFLOW软件非结构网格亚跨声速湍流模拟的验证与确认

计算流体力学(computational fluid dynamics,CFD)数值模拟在航空航天等领域发挥越来越重要的作用,然而CFD数值模拟结果的可信度仍然需要通过不断地验证与确认来提高.本文给出了从制造解精度测试、简单到复杂外形湍流模拟网格收敛性研究等三个方面开展CFD软件验证与确认的方法,并对自主研发的CFD软件平台HyperFLOW在非结构网格上模拟亚跨声速湍流问题的能力进行了验证与确认.首先通过基于Euler方程和标量扩散方程的制造解精度测试,分别验证了HyperFLOW在非结构网格上对Euler方程和黏性项的求解精度,结果表明其能够在任意非结构网格上达到设计的二阶精度. 其次,通过NASATurbulence Modeling Resource中的湍流平板、二维翼型近尾迹流动、二维Bump等几个典型的亚声速湍流算例的网格收敛性研究,量化考察了数值结果的观测精度阶和网格收敛性指数,并与国外知名CFD解算器CFL3D,FUN3D的计算结果进行了对比,验证了HyperFLOW对简单湍流问题的模拟能力,且具有良好的网格收敛性和计算精度(阶). 最后,通过NASA CommonResearchModel标模定升力系数的网格收敛性研究和升阻极曲线预测,验证了软件在复杂外形亚跨声速湍流流动数值模拟中也具有良好的可信度.      

等热流密度加热垂直下降液膜的研究

利用等热流密度加热条件下降膜流动的三维模型方程进行线性稳定性分析和数值模拟。线性稳定性分析表明，模型方程在小到中等Reynolds数下都适用，并且流向不稳定性增长率随着Reynolds数和Marangoni数增加而增加，展向不稳定性增长率则随着Marangoni数增加而增加，随着Reynolds数增加而减小，流向和展向对扰动波数都存在一个不稳定区间。三维数值模拟表明，在等热流密度加热条件下，液膜在随机扰动的情况下最终会形成带孤立波的三维溪流状结构，液膜与气体的换热也因溪流状结构的出现而加强；在随机扰动的基础上引入占优势地位的展向最不稳定扰动会使得换热增强，液膜会提前破裂；在随机扰动的基础上引入占优势地位的流向最不稳定扰动时，液膜的换热会增强，但不会提前破裂；在随机扰动的基础上同时引入占优势地位的流向和展向最不稳定扰动时，换热会加强且液膜会提前破裂。     

一种考虑界面相物性影响的界面单元法

界面是由复杂的界面相简化而成的,界面破坏实际是界面相材料的破坏。数值计算为了方便,如经典模型和内聚力模型等,都把很薄的界面相作无厚度化处理。导致只能考虑界面的面力,而无法考虑界面相内的应力(平行于界面方向的应力)。使界面失效准则先天性地排除了界面相内部应力的影响,从界面相材料失效机理的角度来看这是不够严谨的。本文将界面相材料等效为一种弹性连续体,由界面本构关系推导得到了一种新的界面单元。该单元具有界面参数易确定、对界面相物性可以进行等效描述等优点。通过商用有限元软件ABAQUS和用户子程序UEL实现了数值分析,并与直接物理模型的数值模拟结果进行对比,证明了本方法的简便及准确性。通过对不同界面相厚度结构的进一步分析,探讨了本文方法的可行范围。     

Molecular insight into structures of monocationic and dicationic ionic liquids in mica slits

Structures of ionic liquids (ILs) 1-decyl-3-methylimidazolium bis(trifluoromethanesulfonyl)azanide ([C₁₀mim][TFSA]) and 1-decyl-dimethylimidazolium bis(trifluoromethanesulfonyl)azanide ([C₁₀(mim)₂](TFSA)₂) in different-sized mica slits have been investigated using molecular dynamics simulations. Ion density and angular distributions for monocationic IL [C₁₀mim][TFSA] were analysed to elucidate the IL structures under different surface charges and especially their changes in the direction perpendicular to the surfaces. [C₁₀mim][TFSA] formes in bilayers, compatible with existing models of ILs with long alkyl chains. For dicationic IL [C₁₀(mim)₂](TFSA)₂, cations adjacent to the mica surface tend to stay parallel to the surface with both positively charged rings absorbed. While near the centre of the slit, dications show the weak tendency of orientation distribution, more random than [C₁₀mim]⁺ ions. Structures of [C₁₀(mim)₂](TFSA)₂ cannot be described by bilayer models. Additionally, the in-plane arrangement of [C₁₀mim][TFSA] is more ordered when K⁺ ions completely neutralise the negative charge of the mica surface, and [C₁₀mim]⁺ ions tend to be located in hexagonal mica lattices with two aluminium atoms in replacement of silicon atoms. [TFSA]^? ions are constrained by the neighbouring K⁺ ions absorbed onto mica lattices.